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Abstract 

We explore spin and charge transport phenomena in two dimensional elec- 
tron gas in presence of Rashba spin-orbit coupling connected to two ideal 
Ferromagnetic leads. In particular we show through a combination of analyt- 
ical and numerical calculation that the spin polarization which is transported 
depends on the Magnetization direction of ferromagnet even if the magneti- 
zation of both FM's are parallel. Conductance is also shown to be anisotropic. 
These anisotrpies present in spin and charge transport are a consequence of 
breaking of rotational invariance due to Rashba spin-orbit interaction and 
are present irrespective of the Hamiltonian considered being an effective mass 
Hamiltonian or tight binding model Hamiltonian. 
PACS numbers: 72.25-b,72-25.Dc,72.25Mk,72.25Rb,72.25Hg 

The growing field of spintronics has attracted a lot of interest after the proposal of spin- 
filed effect transistor (SPIN-FET) by Datta and Das 0. The Datta-Das SPIN-FET is a 
hybrid structure of type FM1-2DEG-FM2, where 2DEG is a two-dimensional electron gas 
of a narrow gap semiconductor (InAS) and FM1 and FM2 are injector and detector Fer- 
romagnetic contacts. The working of this device relies on the manipulation of electronic 
spin state in 2DEG with the electric field of an external gate electrode. Essential for this 
mechanism is field dependent spin-orbit coupling, which is relatively large and well estab- 
lished ||. It is now generally accepted that the spin-orbit coupling in narrow-gap 2DEG is 
governed by Rashba Hamiltonian ||. For a 2DEG lying in xy plane(see Fig.l) the Rashba 
spin-orbit interaction has the form, Hr = a(k x a) ■ z,with k being momentum vector, a 



Pauli matrices and z is a unit vector perpendicular to 2DEG plane. The Rashba spin-orbit 
causes spin splitting for k / 0, AE = 2ak, which is linear in momentum. The Rashba 
splitting is due to lack of space inversion symmetry and not due to lack of time reversal 
symmetry. Since Rashba Hamiltonian is time reversal invariant. However the exchange 
splitting in Ferromagnets is due to the breaking of time reversal symmetry. Therfore it is 
natural to expect that spin and charge transport properties of a hybrid structure like SPIN- 
FET, which combines elements with different symmetry properties, may be different than 
the standard mesoscopic structures consisting of elements with same symmetry, for, e.g., all 
metal mesoscopic structures. 

Motivated by this, in this paper we study the spin and charge transport of a FM1-2DEG- 
FM2 system shown in FIG. 1. The question we are addressing is the following: Consider 
the Fig.l, a natural reference frame for the Fig.l is defined by the plane of 2DEG (we 
call it xy plane) and the normal to this plane , i.e., the z axis. The polarization of the 
Ferromagnets FM1 and FM2 are parallel to each other and points in a direction (9, <p) with 
respect to the natural coordinate system i.e. make an angle 9 with z axis and an angle 
with the x axis. Now the question is does the spin polarization which is transported 
from FM1 to FM2 and the charge transport, i.e, conductance depends on (#,</>)? In other 
words if we rotate the polarization vector of Ferromagnets simultaneously with respect to the 
natural coordinate system in such a way that they always remains parallel, does conductance 
and spin polarization which is transported changes? Naively speaking one would expect 
that the conductance and transported spin-polarization should be independent of (9, <fi) 
as long as both the Ferromagnets are parallel. In contrast to naive expectation we show 
through a combination of analytical and numerical calculation that spin polarization which 
is transported and charge conductance are anisotropic and these anisotropics are present 
irrespective of the Hamiltonian considered being an effective mass Hamiltonian Q or tight 
binding Hamiltonian [0. In this sense this is a rather general principal which says the 
polarization of transported electron across a FM/2DEG/FM2 and conductance is anisotropic 
and is a consequence of breaking of rotational invariance due to spin-orbit coupling. This is 
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in contrast to the claim made by Molenkamp et. al. |4| that the effective mass Hamiltonian 
does not have conductance anisotropy while the tight binding model || has due to reduced 
symmetry of lattice. Another important consequence of our study is it points out that spin 
coherence is also anisotropic, i.e. it depends on the chosen basis. 

The Hamiltonian of a 2DEG lying in xy plane (as shown in Fig. 1), in presence of Rashba 
spin-orbit coupling reads: || 

hv 2 

H = --z— + «(axk)-z (1) 
2m 

where a Rashba spin-orbit interaction parameter and a=(a x ,a y ,a z ) is the vector of Pauli 
matrices and z is unit vector along the z axis. We write the above Hamiltonian in the matrix 
form which is more convenient for the study of spin transport 

H=^(B I + B R -a) (2) 

where I is the 2x2 identity matrix, B = — * y and the vector is ~Q^=2a{k y 'St — k x y). 

An appropriate physical quantity to study the spin transport is the Polarization vector 
P=< a > where angular brackets represents the ensemble averaging. With this definition 
one can immediately write down the equation of motion for polarization vector, 

dP d< a > i 

— = : = — — < aH - Ha > 

at at n 

= ^-< ( j(B R - ( T)-(B R - ( T) ( r> (3) 

simplifying above equation using vector identities for triple product and commutation rela- 
tion for Pauli matrices leads to following equation of motion for polarization vector, 

- = B R x<a>. (4) 

The eq. (|J) is well know in the literature and is a fully quantum mechanical and holds eve 
if Br is time dependent. The eq. (||) can be solved analytically when the field Br is a 
constant vector, the most general solution is given as; 
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P(t) = P cos{cu R t) + 2Br(B r ■ P ) sin 2 (^t/2) 

+ (B R x P ) sin(w fl t) (5) 

where Po is the initial polarization imposed by Ferromagnet FM1 (we are interested in 
the case when the Polarization of both the Ferromagnets FM1 and FM2 are parallel and 
equal in magnitude, i.e, Pi = P 2 = Po) , u R = B R /Ti is precession frequency (precession 
angle = to at), B R = 2a ^k\ + ky is the magnitude of Rashba field Br( the direction 
of Br is always perpendicular to the instantaneous wave vector k). During electrons free 
flight the direction and magnitude of Br remains constant hence the solution provided by 
eq.(|5]) is applicable only during the free flight. Since scattering from impurity or boundary 
changes the momentum and hence the filed Br, so the time occurring in eq.(|5|) is free flight 
time. However for the ballistic transport we will use the eq.([|) and take into account the 
boundary scattering later in diffusive approximation as we will see later. Now since we are 
interested in the transport properties when the polarization of both the Ferromagnets are 
parallel to each other but pointing in arbitrary direction (6, <ft) such that P x = P 2 = Po = 
Po(sin # cos </>, sin 6 1 sin 0, cos ^) such that with respect to hence by projecting P(t) on Po we 
get the, 

P(t) • P = |P | 2 cos(uj R t) + (P • B R ) 2 sinV^/2), (6) 

where u R = B R /h = 2akf/h. The eq.(^) is the quantitative measure of spin polarization 
which gets transported through the 2DEG from FM1 to FM2. For a given injection angle 
(3 as shown in Fig. 1, the eq.(f|) simplifies to, 

Pol(0,<P,L,W,u; R t) = P ® ' P ° 

I -to I 

= cos(uj R t) + sin(/3 - 0) 2 sin(#) 2 sm 2 (u R t/2) (7) 

In the above equation t is the time electron takes to reach the output terminal and it is clear 
from eq.([F[) that value of transported polarization ,i.e, Pol(9, (ft, L, W, ou R t), lies between +1 
and -1. Since the electron are injected over — 7r/2 < (3 < n/2, we need to make an average 
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over all possible values of injection angle (3. However depending upon injection angle (3 elec- 
tron reaches the boundary without scattering(dashed trajectory in Fig. 1) or with scattering 
(solid trajectory in Fig.l) from the boundaries. Hence we need to calculate t accordingly 
for different values of (3. Therefore we divide the the integration over j3 in three regimes, 
namely (a) -vr/2 < (3 < -tan- 1 (W/2L), (b) - t&ir 1 (W/2L) < (3 < - taxr 1 (W / 2L) and 
(c) tan~ 1 (H / /2L) < (3 < it/ 2 where (a) and (c) corresponds to the trajectories which suf- 
fers scattering from boundary while trajectories in regime (b) propagates ballistically. For 
injection angle (3 in the regime (b) electron reaches the output terminal(FM2) ballistically 
therfore the time to reach the output terminal is t = L/ cos(/5)(see Fig.l dashed line), while 
for regime (b) and (c) electron scatters from the boundary at least once before reaching the 
out put terminal(FM2), hence for these values of (3 we assume that the electrons diffuse 
along the channel with a mean free path W/2 sin(/?) (later in our exact numerical simulation 
we will see that this approximation is quite reasonable) so the time to reach the boundary 
is given as t = (2L 2 sm(/3))/(vfW). Using these vale for t we get 



u R t 



2ak f L 

v f cos(/3) — cos(/3) ^ ^ L V2L^' V 2L 



i±- = Me<± p e ^tan-^f), tan-^f )} 



(8) 

v f W = W y ^ l-"-2' ua,±1 



2afcfL2sin(/3) = 4 "° L " 2sin(/3) (3 e {±f , itan-^f )} 



where a = akf/Ef is dimensionless Rashba parameter(.Ey is Fermi energy), L = L/Xf and 
W = W/Xf is length and width of channel in units of Fermi wavelength. Substituting 
these values of u^t in eq. (0) and performing the integration over f3 we obtain polarization as 
function of 9, 0, L, W, a. Eq.(0) together with eq.(||) can be used to calculate the transported 
polarization for any given direction (6,<p), however for clarity and simplicity we present 
results for three specific cases corresponding to different values of 9 and 0, namely, (i)#=0, 
is variable, i.e., polarization of FM1 nad FM2 is rotated in xy plane (the plane formed 
by 2DEG) (ii) 0=0, 9 is variable corresponding to the rotation in xz plane (iii)0=7r/2, 9 
is variable, corresponding to the rotation in yz plane. For these three different cases the 
transported polarization given by eq.(|7]) is shown in Fig. 2 as a function of angle , where the 
other parameters are L = W = 50/(27r) and a = 0.06. It is clearly seen from Fig. 2 that 
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polarization which is transported is anisotropic, it is a consequence of spin-orbit coupling 
which breaks the rotational symmetry. The amplitude of oscillation tells about the spin 
coherence and since this is different for all the three cases, signifying that the spin coherence 
is also affected anisotropically. Infact it is seen from fig. 2 that amplitude of oscillation is 
larger for the case (i) and (ii), when the polarization vector of Ferromagnets lies in yz or 
xy, compared to the case (iii). The absolute magnitude of oscillation is always smaller than 
one implying even in ballistic transport spin dephasing takes place due to the boundary 
scattering . 

To further strengthen our results we performed numerical simulation on a tight binding 
square lattice of lattice spacing a with N x sites along x axis and N y sites along y axis. For 
tight binding Hamiltonian the Rashba spin-orbit coupling is given by X so = a /2a = atkfa/2. 
We fix t=l (hopping) and kfa=l (ballistic case) for numerical simulation in tight binding 
model. Once t and kfa are fixed the other parameters for tight binding model which would 
corresponds to the parameters of Fig. 1 are given as, N x = 2ttL = 50,N y = 2nW = 50 and 
X so = atkfa/2 = 0.03. With these set of parameters we calculate spin resolved conductance 
for a given polarization direction (8, 0) of Ferromagnets, within Landauer-Biittiker formalism 
Using the spin resolved conductance we define polarization as 

P = ^^ L - 0) 

<^sc + <**sf 

where G sc and G s f are spin-conserved and spin flip conductance respectively. The quantity 
P in eq.(||) corresponds to the quantity given in eq.(0) and also lies between +1 and - 
1. This is plotted in Fig. 4 , we see that the agreement between Fig. 3, i.e, analytical 
calculation, and Fig. 4 is quite good. The slight quantitative mismatch is due to the fact 
that numerical simulation was done for hard wall confining potential in y direction which 
leads to specular reflection, while in analytical calculation scattering from the boundary was 
treated as diffusive. Therefore it is clear that the anisotropy in spin transport is present in 
continuum model (effective mass Hamiltonian) as well as in tight binding model and is not 
an effect of reduced symmetry of tight binding model, rather it is a consequence of breaking 
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of rotational invariance due to spin orbit coupling. 

Now since conductance of FM/2DEG/FM, depends on the polarization of electrons reach- 
ing the output terminal, hence it is expected that conductance should also be anisotropic. 
This is clearly visible in Fig. 5 where we have plotted total conductance, i.e., G = G sc + G s f 
corresponding to the Fig. 4, as function of polarization angle. It should be noted that the 
conductance is symmetric with respect to angle 9 or <fi which is consistent with Buttiker 
symmetry relation for charge transport ||. It is important to point out that in recent lit- 
erature H an erroneous result was reported, where it was claimed that conductance of a 
FM/2DEG interface changes on flipping the magnetization of FM which is incorrect. 

The results presented above were in ballistic regime. To verify that these results survives 
in diffusive case we show polarization and conductance in Fig. 5 and Fig. 6 respectively for 
diffusive case. We have taken Anderson model for disorder with width 3|t|, corresponding to 
a mean free path of I = 10a. The other parameters are same as those for Fig. 3 and Fig. 4. 
It is clearly seen that the anisotropy survives even in diffusive case. This only strengthen our 
previous assertion that spin coherence is anisotropic. Also it is instructive to compare Fig. 
3 for ballistic transport and Fig. 5 for diffusive transport. It is seen that the polarization 
which is transported is not affected much by the presence of disorder which is consistent 
with the Rashba spin-orbit interaction which is independent of disorder strength. However 
the magnitude of charge conductance is reduced drastically as seen from Fig. 4 and Fig. 6, 
though the qualitative behavior as function of angle remains unchanged. This clearly demon- 
strates that the conductance anisotropy exist and is consistent with the Buttiker symmetry 
relation. One important thing to be noticed is the amplitude of oscillation for ballistic case 
as well for diffusive case for both polarization and conduct anceremains almost unchanged 
since the Rashba coupling was kept fixed for all the figures. This clearly demonstrates that 
the anisotropy is a consequence of spin-orbit interaction and is not affected by disorder. 

I thank P. Bruno and G. Bouzerar for helpful discussion. 
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FIGURE CAPTINOS 



lFig.l A 2DEG connected to two ideal Ferromagnetic leads. 

2Fig2. Polarization as a function of angle calculated using eq.(^) and eq.(^j). Where L = 
W = 50/2vr,a = 0.06. 

3Fig.3 Results of numerical simulation for polarization for ballistic system. The numerical 
simulation were performed on a 50 x 50 lattice within tight binding model. The tight 
binding Rashba parameter is given by X so = atkfa/2 = 0.03, FM exchange splitting 
is A/Ef = 0.5 and kfd=l. These parameters were chosen in such a way that they 
correspond to the parameters of Fig. 1, as explained in text. 

4Fig.4 The conductance as a function of angle. The parameters are same as in Fig.3 

5 Fig. 5 Polarization as a function of angle for diffusive case. Here kfl—10, where / is mean 

free path. Configuration averaging was performed over 15 different configuration. The 
other parameters are same as in Fig. 3 

6 Fig. 6 Conductance as function of angle corresponding to the Fig. 5. Here kfl=10, 

where / is mean free path. Configuration averaging was performed over 15 different 
configuration. The other parameters are same as in Fig. 3 
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